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Abstract 



> 

Recent experimental results and developments in the theoretical treatment of neutrino-nucleus inter- 
5^ actions in the energy range of 1-10 GeV are discussed. Difficulties in extracting neutrino-nucleon cross 
sections from neutrino-nucleus scattering data are explained and significance of understanding nuclear 
effects for neutrino oscillation experiments is stressed. Detailed discussions of the status of two-body 
current contribution in the kinematic region dominated by quasi-elastic scattering and specific features 
of partonic nuclear effects in weak DIS scattering are presented. 
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0.1 Introduction 

Recent interest in neutrino interactions in the few GeV energy region comes from neutrino osciUation 
experiments and their need to reduce systematic errors. Neutrino fluxes used in contemporary long 
and short basehne experiments (K2K, T2K, MINOS, NOvA, MiniBooNE) are peaked in the 1 - 5 GeV 
energy domain and during the last ~ 10 years there has been considerable theoretical and experimental 
activity in the investigation of neutrino cross sections in this domain with reference [1] being a good 
summary of the lower-energy situation. Several new cross section measurements have been performed 
by neutrino oscillation collaborations and there are two dedicated cross section experiments (SciBooNE 
and MlNERvA) which have been launched at Fermilab. 

Even with this degree of activity, the precision with which the basic neutrino- nucZeon cross sections 
are known is still not better than 20 — 30%. There are two main reasons for this: the poor knowledge 
of neutrino fluxes and the fact that all the recent cross section measurements have been performed 
on nuclear targets. It is important to recall that what current neutrino experiments are measuring 
are events that are a convolution of energy-dependent neutrino flux (8) energy-dependent cross section 
(8) energy-dependent nuclear effects. The experiments have, for example, then measured an effective 
neutrino-carbon cross sections and to extract a neutrino-nucleon cross sections from these measure- 
ments requires separation of nuclear physics effects that can be done with only limited precision. For 
many oscillation experiments, using the same nuclear targets for their near and far detectors is a good 
start. However, even with the same nuclear target near-and-far, that there are different near and far 
neutrino energy spectra due to oscillations means there is a different convolution of cross section 
nuclear effects near and far and there is no automatic cancellation between the near-and-far detectors. 
For a thorough comparison of measured neutrino-nucleon cross sections with theoretical models, these 
convoluted effects have to be understood. 

Some of the new cross section measurements raised doubts in the areas which seemed to be well 
understood. The list of new puzzles is quite long and seems to be expanding. What is the value of 
the quasielastic axial mass? How large is the two-body current contribution that can mimic genuine 
quasielastic interactions? How large is CC (charged current) coherent pion production at a few GeV 
neutrino energies? What is behind the large discrepancy between MiniBooNE pion production mea- 
surements and theoretical model predictions? It can be seen as a paradox that the more than 30-year 
old ANL and BNL low statistics deuterium pion production data, with its minimal nuclear corrections, 
is still used as the best source of information about the nucleon-A transition matrix element. 

Analysis of neutrino scattering data is certainly more complicated than the analysis of electron 
scattering data. In the electron case one knows exactly the initial electron energy and so also the 
values of energy- and momentum-transfer. It is then possible to explicitly study separate interesting 
kinematical regions like QE (quasielastic) peak or the A peak. Neutrino scattering data is always flux 
(often wide band!) integrated. Interacting neutrino energy must be evaluated based on kinematics of 
particles in the final state taking into account detector acceptance and measurement accuracy. 

For neutrino- nwc/eon interactions one can distinguish: Charged Current Quasielastic (CCQE), 
Neutral Current elastic (NCEl), Resonance production (RES) and more inelastic reactions up to the 
deep-inelastic (a rather misleading "DIS" term is often used to describe all the interactions which are 
neither CCQE/NCEl nor RES) domain. Quite different theoretical tools are used to model each of 
them. The simplest neutrino hadronic reaction is the charge current quasielastic (CCQE) interaction: 
I'i + n ^ i~ + p with two particles: charged lepton and proton in the final state. One would like to 
extend this definition to the neutrino-nucleus interaction occurring on bound neutrons. The obvious 
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question arises: what is the experimental signature of CCQE on a nuclear target? The ejected pro- 
ton is not necessarily seen in a detector because quite often its momentum is below the acceptance 
threshold. However, events with a single reconstructed charged lepton track can result from a variety 
of initial interactions eg. from a two body current interaction or from real pion production and its 
subsequent absorption. Similar problems arise in other type of interactions. It is becoming clear that 
interpretation of neutrino-nucleus interaction must rely on a careful data/Monte Carlo (MC) compar- 
ison done with reliable MC neutrino event generators. This is why we decided to include in the review 
some information about development of MC event generators. 

From the experimental point of view it is natural to speak about events with no pions in the final 
state, with only one pion etc. In fact, in several recent experimental measurements that investigated 
quantities defined in this way, the dependence on assumptions of Monte Carlo event generators were 
minimal. To compare with the experimental data given in this format one must add contributions from 
various dynamical mechanisms and also to model FSI effects. Several ingredients of the theoretical 
models are verified simultaneously. It is clear that in order to validate a model one needs many samples 
of precise neutrino-nucleus scattering measurements on variety of nuclear targets with various neutrino 
fiuxes. 

Our review is organized as follows, we review recent inclusive measurements in the lower E region 
and then concentrate on exclusive states in increasing W, the mass of the hadronic system. Due to the 
limited length of this review, we do have to limit our coverage to only the most recent developments. 

0.2 Neutrino Charged Current and Neutral Current Inclusive Re- 
actions 

0.2.1 Recent measurements 

There are four recent CC inclusive neutrino and antineutrino cross sections measurements in the 
E,^ < 10 GeV energy region [5], see Fig. [l| We notice a mild tension between SciBooNE and T2K 
measurements. In the following sections QE, RES and DIS contributions will be discussed separately. 



0.2.2 Theory. General formulae: outgoing lepton differential cross sections 

In this paper, we will discuss the neutrino CC or NC (neutral current) inclusive reaction: 

ueik) + Az ^ £{k') + X. (1) 

The generalization of the expressions to antineutrino induced reactions is straightforward. In the 
equation above, the outgoing lepton could be either a negatively charged lepton, £~ , of fiavor £ or a 
neutrino for CC or NC processes, respectively. 

The double differential cross section, with respect to the outgoing lepton kinematical variables, for 
the process of Eq. ^ is given in the Laboratory (LAB) frame by 

d'a^.i _ \k'\ (2) 



d^{k')dE[ |A?|47r2r/ 

with k and k' the LAB lepton momenta, E[ = {k'"^ + rnj)^^"^ and me the energy and the mass of the 
outgoing lepton, Gp = 1.1664 x 10^^^ MeV"^, the Fermi constant and L and W the leptonic and 
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Figure 1: Breakdown of the inclusive CC muon neutrino cross section on free isoscalar target to QE, RES and DIS 
contributions, as viewed by NuWro MC event generator. 

hadronic tensors, respectively. Besides, rj takes the values 1 or 4 for CC or NC processes, respectively. 
The leptonic tensor is given by (in this convention, eoi23 = +1 and the metric is g^'^ = (+, — , — , — )): 

L^^J = L^^ + iL^ij = k^ka + k^k^ — g^ak • k + it^aapk k^ (3) 

The hadronic tensor includes a collection of non-leptonic vertices and corresponds to the charged or 
neutral electroweak transitions of the target nucleon or nucleus, i, to all possible final states. It is thus 
given by 

= ^E(2^)''^'(^/-^-9)(/bcc+,nc(0)N>(/|jTc+,nc(0)N>* (4) 

with the four-momentum of the initial target, = the target mass square, P'j the total four 
momentum of the hadronic state / and q = k — k' the four momentum transferred to the hadronic 
system. The bar over the sum denotes the average over initial spins. 

The hadronic tensor is completely determined by six independent, Lorentz scalar and real, structure 
functions Wi{q^ , Q ■ P), 

Taking g in the z direction and = (Mj,0), it is straightforward to find the six structure functions 
in terms of the W^^,W^^ = ,W^^ ,W^y and W^'"^ components of the hadronic tensor. After 
contracting with the leptonic tensor, one obtains that for massless leptons only three of them are 
relevant, namely 



dxdy r]TT 



{ (l - y - ^) F,^ + xy'Fi' + yil - y/2)xF.^^ (6) 
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with E,y the incoming neutrino energy, M the nucleon mass, x = — q^/2Mq^, y = q^/E^, while the 
nuclear structure functions F^2,3 are given by, = 2MMiWi, = 2{q ■ P)W2 and F^/M = 
—2{q - P)W'i / Mi. The cross section for the CC antineutrino induced nuclear reaction is easily obtained 
by i) changing the sign of the parity-violating term, proportional to F3, in the differential cross section 
(this is because L^^} = L^/l.), Eq. (6), and ii) using — Jcc+ the definition/computation of the 
hadron tensor in Eq. Q. In the case of antineutrino NC driven processes, it is only needed to flip 
the sign of the term proportional to -F3 in the differential cross section, since the hadron NC is not 
affected. 

The hadronic tensor is determined by the W or Z gauge boson selfenergy, Il'^^^(q), in the nuclear 
medium. To evaluate this latter object requires a theoretical scheme, where the relevant degrees of 
freedom and nuclear effects could be taken into account. 

In the next two sections we will discuss CCQE and pion production reaction. The general formalism 
described above will be used in the section devoted to DIS. 



0.3 Charged Current Quasielastic 

As discussed in the Introduction, we define CCQE as either the reaction on a free nucleon or on a quasi 
free nucleon inside a nucleus yielding a muon and nucleon. In the case of neutrino nucleus scattering 
we also use the term CCQE-like reaction defined as one in which there are no pions in the final state. 
It then includes events with real pion production followed by absorption. Such a definition may seem 
awkward but as will be seen, it is close to what was experimentally measured by the MiniBooNE 
collaboration. 

A theoretical description of the free nucleon target CCQE reaction is based on the conserved 
vector current (CVC) and the partially conserved axial current (PCAC) hypotheses. The only un- 
known quantity is the nucleon axial form-factor Ga(Q^) for which one typically assumes a dipole form 
Ga(0)(1 + t7t)~^ with one free parameter, the axial mass Ma- The non-dipole axial form factor was 
investigated e.g. in [3j. 

In the past, several measurements of M^ were performed on a deuterium target for which most 
of nuclear physics complications are minimal and it seemed that the results converged to a value 
of the order of 1.03 GeV There is an additional argument in favor of a similar value of Ma_ 
coming from the weak pion-production at low Q^. PCAC based evaluation gives an axial mass value 
of 1.077 lb 0.039 GeV [5]. On the other hand, all of the more recent high statistics measurements 
of M4, with the exception of the NOMAD higher-energy experiment, reported larger values: K2K 
(oxygen, > 0.2 GeV^) ^ 1.2 ± 0.12 [6]; K2K (carbon, > 0.2 GeV^) ^ 1.14 ± 0.11 [7]; MINOS 
(iron, Q2 > GeV^) ^ 1.19 ± 0.17; MINOS (iron, > o.3 GeV^) ^ 1.26 ± 0.17 [8]; MiniBooNE 
(carbon, > q GeV^) 1.35 ±0.17 [9j; MiniBooNE (carbon, > 0.25 GeV^) 1.27 ±0.14 (for 
completness: NOMAD (carbon, > q GeV^) 1.07 ±0.07 [10]). 

The difference between MiniBooNE and NOMAD measurements could come from different def- 
initions of the CCQE signal. In the case of MiniBooNE a sample of 2-subevents (Cherenkov light 
from muon and from decay electron) is analyzed and ejected protons are not detected. In the case 
of NOMAD 1-track (muon) and 2-tracks (muon and proton) samples of events are analyzed simu- 
lateuosly. With a suitable chosen value of the formation zone parameter tq values of M^ extracted 
separately from both data samples are approximately the same, see Table 9 in [10] . We note that the 
procedures in which the formation zone concept is applied to nucleons that already exist may seem 
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little controversial. We would like to mention also the CCQE data not yet published in peer review 
journals. MINOS tried to evaluate better the pion production background [11]. A function of 
which corrects Monte Carlo (NEUGEN) RES predictions was proposed. The shape of the curve is 
similar to MiniBooNE's DATA/MC correction function (see below) but in the case of MiniBooNE for 
Q2 > 0.1 GeV^ the correction factor is > 1. The new MINOS best fit value of Ma is 1.16 GeV and the 
error was reduced by a factor of 3 with respect to [S]. SciBooNE showed partial results of the CCQE 
analysis [12]. Results are given in terms of fits for CCQE cross-section DATA/MC multiplicative 
factors aj (j label true neutrino energy bins) and a scaling factor Fjv- The obtained best fit values in 
the neutrino energy region Ey G (0.6, 1.6) GeV are between 1.00 and 1.09 which with F]\[ = 1.02 and 
the value of the axial mass used in the NEUT Monte Carlo generator (1.2 GeV) should translate to 
the axial mass value Ma ~ 1.25 — 1.3 GeV. In the SciBooNE analysis there are some instabilities in 
the wider region of (see Fig. 11.2 in [13]). A use of the universal background scaling factor cibcg 
for three different event samples is perhaps not sufficient (its best fit value is as large as 1.37 GeV). 

An important antineutrino CCQE measurement was reported by MiniBooNE [TJj. The DATA/MC 
average cross-section ratio was reported to be 1.21 it 0.12 which is a surprising result because in the 
NUANCE carbon CCQE computations the Ma value was set to be 1.35 GeV. In the experimental 
analysis, it was important to evaluate correctly neutrino contamination in the anti-neutrino fiux. Three 
independent measurements indicate that the fiux in the antineutrino beam should be scaled down 
by a factor of ~ 0.8 with an obvious important impact on the final results. 

The most recent MINERvA preliminary results for CCQE antineutrino reaction are still subject 
to large fiux normalization uncertainties but they seem to be consistent with Ma = 0.99 GeV |15j . 

0.3.1 MiniBooNE data 

In recent discussions of the CCQE, MiniBooNE measurement plays a special role. For the first time 
the data was presented in the form of double differential cross section in muon scattering angle and 
kinetic energy. Such data is the actual observable for the MiniBooNE experiment and more complete 
than a distribution of events in which is calculated assuming an obviously incorrect nuclear model 
(the nucleon is assumed to be at rest). The signal events form a subset of events with no pions in the 
final state. MiniBooNE subtracted as a background, events with real pion production and subsequent 
absorption and also a contribution from pionless A decays implemented in the NUANCE MC [16] 
as constant fractions of and A"*" decays, following the approach of Ref. |I7j . The background 
estimate, based on MC predictions, was later corrected by a dependent function, which accounts 
for a data/MC discrepancy in the sample of events containing one 7r+ in the final state. The shape of 
the correction function is not well understood [T8] but it has an important impact on the extracted 
value of Ma- The function quantifies a lack of understanding of processes like pion absorption and 
can have a significant eff'ect on the understanding of both samples of events. 

MiniBooNE also provided data for the CCQE signal plus background together as the measurement 
of the cross section of the process in which there are no pions in the final state, the observable which 
is maximally independent of MC assumptions. 

0.3.2 Theoretical approaches to CCQE - generalities 

Several approaches have been followed/derived to compute the relevant gauge boson absorption modes 
(self-energy) to describe the CCQE process. For moderate and intermediate neutrino energies, in the 
few GeV region, the most relevant ones are: the absorption by one nucleon, or a pair of nucleons or 
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even three nucleon mechanisms, real and virtual meson (vr, p, ■ ■ ■) production, excitation of A of higher 
resonance degrees of freedom, etc. (for example, some absorption modes are depicted in Fig. [2] for the 
case of neutrino CC processes). A review of theoretical model results can be found in [T^]. Almost all 
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Figure 2: Diagrammatic representation of some diagrams contributing to the H^^— selfenergy. 

approaches, used at intermediate neutrino energies, deal with hadron, instead of quarks and gluons, 
degrees of freedom. In addition they consider several nuclear effects such as RPA or Short Range 
Correlations (SRC). The free space couplings between hadrons and/or the weak W and Z bosons are 
parametrized in terms of form factors, which are fitted to the available data on electroweak scattering 
off free nucleons. In the few GeV energy region, theoretical models rely on the impulse approximation 
(lA) and neutrino-nucleus CCQE interactions are viewed as a two step process: primary interaction 
and Final State Interactions (FSI), propagation of resulting hadrons through the nucleus. The validity 
of the I A is usually related to typical values of the momentum transfer q. Experience from the electron 
scattering tells us that for q > 300 — 500 MeV/c lA based models are able to reproduce the data well. 
Thus, the expectations is that for a few GeV neutrino interactions lA is an acceptable approach 
and if necessary simpler nuclear models computations can be supplemented with RPA corrections for 
lower momentun transfers (see below). In the neutrino nucleus cross section measurements a goal is 
to learn about neutrino free nucleon target scattering parameters (an obvious exception is coherent 
pion production). Effective parameters like sometimes discussed quasi elastic axial mass are of 

little use as their values can depend on the neutrino flux, target and perhaps also on the detection 
technique/acceptance. 

The definition of neutrino-nucleus CCQE scattering can be made more rigorous in the language 
of many body field theory. CCQE process originates from a first step mechanism where the gauge 
boson is being absorbed by just one nucleon. This corresponds to the first of the selfenergy diagrams 
depicted in Fig. [2] (contribution (a)). This contribution, that from now on we will call genuine QE, 
has been computed within different theoretical models and used to predict the corresponding outgoing 
lepton differential cross section. 
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The simplest model, commonly used in Monte Carlo event generators, is the relativistic Fermi 
gas (RFG) model proposed by Smith and Moniz more than 35 years ago |20] corresponding to only 
one many body Feynman diagram. The model combines the bare nucleon physics with a model to 
account for Fermi motion and nucleon binding within the specific nucleus. The model can be made 
more realistic in many way^ to achieve better agreement with a broad range of electron scattering 
data. For example, the inclusion of a realistic joint distribution of target nucleon momenta and 
binding energies based on short range correlation effects leads to the spectral function (SF) approach. 
Spectral functions for nuclei, ranging from carbon (A = 12) to iron (A = 56) have been modeled using 
the Local Density Approximation (LDA) [21], in which the experimental information obtained from 
nucleon knock-out measurements is combined with the results of theoretical calculations in nuclear 
matter at different densities, and they have been extensively validated with electron scattering data. 
Calculations by Benhar et al., f22] and Ankowski et al.. [23] show that the SF effects moderately modify 
the muon neutrino differential cross sections, and they lead to reductions of the order of 15% in the 
total cross sections. This is corroborated by the results obtained within the semi-phenomenological 
model (density dependent mean-field potential in which the nucleons are bound) |24] employed within 
the GiBUU model to account for these effects. 

Inclusion of nucleon-nucleon long-range correlations leads to RPA (Random Phase Approximation) 
which improves predictions at lower momentum transfers (and also low Q"^). RPA corrections have been 
discussed by many authors in the past and recently included in computations of three groups (IFIC, 
Lyon and Aligarhj^ in Refs. [251 [26], 123 EH], and [29j respectively. When the electroweak interactions 
take place in nuclei, the strengths of electroweak couplings may change from their free nucleon values 
due to the presence of strongly interacting nucleons. Indeed, since the nuclear experiments on /3 decay 
in the early 1970s |30j . the quenching of axial current is a well-established phenomenon. The RPA 
re-summation accounts for the medium polarization effects in the Iplh contribution (Fig. [2][^a)) to the 
W and Z selfenergy by substituting it by a collective response as shown diagrammatically in the top 
left panel of Fig. [3] Evaluating these effects, requires an in-medium baryon-baryon effective force, 
which in both sets (IFIC and Lyon) of calculations was successfully used/tested in previous works on 
inclusive nuclear electron scattering. RPA effects are important as can be appreciated in the top right 
panel of Fig.|3j In this plot, we show results from both IFIC and Lyon models, presented in Refs. [31j 
and [.32j, respectively for the CC quasielastic v^—^'^C double differential cross sections convoluted 
with the MiniBooNE flux [33]. There, we also see that predictions of both groups for these genuine 
QE contribution, with and without RPA effects, turn out to be in a quite good agreement. Finally, 
it is important to stress also that RPA corrections strongly decrease as the neutrino energy increases, 
while its effects should account for a low deficit of CCQE events reported by several experimental 
groups (see bottom panels of Fig. [s]). Continuum RPA (CRPA) computations for neutrino scattering 
were performed by the Ghent group |34j . 

Other theoretical developments 

In [65 \ I36 t [57| the bound-state wave functions are described as self-consistent Dirac-Hartree solutions, 
derived within a relativistic mean field approach by using a Lagrangian containing a and uj mesons [38]. 

^When the axial mass and electromagnetic form factors are kept unchanged, the inclusion of more sophisticated 
nuclear effects makes the cross section generally smaller with respect to the RFG (relativistic Fermi gas) model. 

^The Aligarh group uses a similar approach to that derived in [25], but with some simplifications that though well 
suited to study the related process of muon capture in nuclei, might not be totally appropriate for the case of larger 
energies and momenta being transferred to the nucleus (see the discussion in [25) ) . 
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OiilyQE: R PA effects Oi.i.v qE: rpa rfferi. 
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Figure 3: Top Left: Set of irreducible diagrams responsible for the polarization (RPA) effects in the Iplh contribution 
to the W or Z self-energies. Top Right: MiniBooNE flux-averaged CC quasielastic Vy, — ^'^C double differential cross 
section per neutron for 0.8 < cos 6^ < 0.9 as a function of the muon kinetic energy. Bottom: Different theoretical 
predictions for muon neutrino CCQE total cross section off ^'^C, as a function of the neutrino energy (left) and (right), 
obtained from the relativistic model of Ref. [25]. In all cases Ma ~ 1.05 GeV. 



This scheme also accounts for some SF effects. Moreover, these models also incorporate the FSI 
between the ejected nucleon and the residual nucleus. The final nucleon is described either, as a 
scattering solution of the Dirac equation [361 in the presence of the same relativistic nuclear mean 
field potential applied to the initial nucleon, or adopting a relativistic multiple-scattering Glauber 
approach [35] . 

The relativistic Green's function model [39j would be also appropriate to account for FSI effects 
between the ejected nucleon and the residual nucleus for the inclusive scattering, where only the 
outgoing lepton is detected. There, final-state channels are included, and the flux lost in each channel 
is recovered in the other channels just by the imaginary part of an empirical optical potential and the 
total flux is thus conserved. 

Another interesting approach starts with a phenomenological model for the neutrino interactions 
with nuclei that is based on the superscaling behavior of electron scattering data. Analysis of inclusive 
(e, e') data have demonstrated that for momentum transfers q >~ 500 MeV/c at energy transfers 
below the QE peak superscaling is fulfilled rather well |40] . The general procedure consist on dividing 
the experimental (e, e') cross section by an appropriate single-nucleon cross section to obtain the 
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experimental scaling function, which is then plotted as a function of a certain scaling variable for 
several kinematics and for several nuclei. If the results do not depend on the momentum transfer q, 
then scaling of the first kind occurs, if there is no dependence on the nuclear species, one has scaling 
of the second kind. The simultaneous occurrence of scaling of both kinds is called superscaling. The 
superscaling property is exact in the RFG models, and it has been tested in more realistic models of the 
(e, e') reaction. The Super-Scaling approach (SuSA) is based on the assumed universality of the scaling 
function for electromagnetic and weak interactions [JT]. The scaling function thus determined from 
(e, e') data is then directly taken over to neutrino interactions [411142] . There are no RPA correlations 
or SF corrections explicitly taken into account, but they may be contained in the scaling function. 
Nevertheless, such approach is far from being microscopic. Moreover, it is difficult to estimate its 
theoretical uncertainties, as for example to what extent the quenching of the axial current, that is due 
to RPA corrections, is accounted for by means of scaling functions determined in (e, e') experiments, 
which are driven by the vector current. 

Theretical models versus MiniBooNE 2D data 

The MiniBooNE data [9] have been quite surprising. Firstly, the absolute values of the cross section 
are too large as compared to the consensus of theoretical models \19\ H3] . Actually, the cross section 
per nucleon on ^^C is clearly larger than for free nucleons. Secondly, their fit to the shape (excluding 
normalization) of the distribution done within the RFG model leads to the axial mass. Ma = 
1.35 lb 0.17 GeV, much larger than the previous world average (~ 1.03 GeV) [5l [10]. Similar results 
have been later obtained analyzing MiniBooNE data with more sophisticated treatments of the nuclear 
effects that work well in the study of electron scattering. For instance, Refs. |44| H5] using the impulse 
approximation with state of the art spectral functions for the nucleons fail to reproduce data with 
standard values of M^. Large axial mass values have also been obtained in ref. [36] where the 2D 
differential cross section was analyzed for the first time using RFG model and spectral function. Similar 
results were obtained in Ref. [47], where the data have been analyzed in a relativistic distorted- wave 
impulse approximation supplemented with a RFG model. 

0.3.3 Multinucleon mechanisms 

A plausible solution to the large axial mass puzzle was firstly pointed out by M. Martinj^et al. [271 EH], 
and later corroborated by the IFIC group [311 B9] . In the MiniBooNE measurement of Ref. [5] , QE 
is related to processes in which only a muon is detected in the final state. As was already discussed 
above, besides genuine QE events, this definition includes multinucleon processes (Fig. gef), where 
the gauge boson is being absorbed by two or more nucleons, and others like real pion production 
followed by absorption (Fig. ^c) and (d)). The MiniBooNE analysis of the data attempts to correct 
(through a Monte Carlo estimate) for some of these latter effects, such as real pion production that 
escapes detection through reabsorbtion in the nucleus leading to multinucleon emission. But, it seems 
clear that to describe the data of Ref. [9], it is necessary to consider, at least, the sum of the selfenergy 
diagrams depicted in Figs. and (e). Those correspond to the genuine QE (absorption by just 

''The papers of Martini et al are based on the older investigation by Marteau et al |3S]. The relavant features of the 
model were known already at the end of 1990s and at that time the goal was to understand better SupeKamiokande 
atmospheric neutrino oscillation signal. 

''Note that the intermediate pion in this term is virtual and it is part of the AA'^ — >■ A'^A'^ interaction inside of the 
nucleus. Indeed, one should consider a full interaction model for the in medium baryon-baryon interaction. 
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one nucleon), and tlie multinucleon contributions, respectively. The sum of these two contributions 
contribute to the CCQE-like cross sectiorj^ 

The inclusion of the 2p2h contributions enables [Sl\ [52] the double differential cross section 
cos 9^ and the integrated flux unfolded cross sectiorj^ measured by MiniBooNE, to be de- 
scribed with values of M4 (nucleon axial mass) around 1.03 it 0.02 GeV 0[TU]. This is re-assuring 
from the theoretical point of view and more satisfactory than the situation envisaged by some other 
works that described the MiniBooNE data in terms of a larger value of Ma of around 1.3-1.4 GeV, 
as mentioned above. 



Similarites and differences between multinucleon ejection models 
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Figure 4: MiniBooNE flux-averaged CC quasielastic u^ — ^'^C double differential cross section per neutron for 0.8 < 
cos 6^ < 0.9, as a function of the muon kinetic energy. Experimental data from Ref. are multiplied by 0.9. In all the 
cases Ma ^ 1.05 GeV. 



As shown in the top panel of Fig. [3| the IFIC group predictions [311, H9] for QE cross-sections agree 
quite well with those obtained in Refs. [271 EEl |32] (Lyon group). However, both above presented 
approaches considerably differ (about a factor of two) in their estimation of the size of the multinucleon 
effects, as can be appreciated in Fig.|4j IFIC predictions, when the 2p2h contribution is included, favor 
a global normalization scale of about 0.9 (see [31j). This is consistent with the MiniBooNE estimate 
of a total normalization error as 10.7%. The IFIC evaluation in |49[ I31j. of multinucleon emission 
contributions to the cross section is fully microscopical and it contains terms, which were either not 
considered or only approximately taken into account in [271EH1E2]- Indeed, the results of these latter 
works rely on some computation of the 2p2h mechanisms for the (e, e') inclusive reaction ([50]), which 
results are simply used for neutrino induced processes without modification. Thus, it is clear that 
these latter calculations do not contain any information on axial or axial-vector contribution^ For 



^Also for simplicity, we will often refer to the multinucleon mechanism contributions, though they include effects 
beyond gauge boson absorption by a nucleon pair, as 2p2h (two particle-hole) effects. 

^We should warn the reader here, because of the multinucleon mechanism effects, the algorithm used to reconstruct 
the neutrino energy is not adequate when dealing with quasielastic-like events, a distortion of the total flux unfolded 



0.3.5 



cross section shape could be produced. We will address this point in Subsect. 

^The evaluation of the nuclear response induced by these 2p2h mechanisms carried out in Ref. [27j is approximated, 
as acknowledge there. Only, the contributions in [27] that can be cast as a A— selfenergy diagram should be quite similar 
to those derived in [49] by the IFIC group, since in both cases the results of Ref. |17] for the A— selfenergy are used. 
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antineutrinos the IFIC model predicts, contrary to the results of the Lyon group, also a sizeable effect 
of 2p2h excitations. 

Another microscopic approach to 2p2h excitations was proposed by Amaro et al. These authors 
have used the empirical (e, e') SuSA scaling function to describe the CCQE MinibooNE data, includ- 
ing some 2p2h contributions due to MEC (meson exchange currents) |51l E2] ■ The approach, used in 
these latter works, to evaluate the 2p2h effects, though fully relativistic, does not contain the axial 
contributions. The authors of |51l [52] also find an increase of the inclusive cross section for neutri- 
nos; at forward muon angles the calculations come close to the data, but the MEC contributions die 
out fast with increasing angle so that the cross section is significantly underestimated at backward 
angles. As a consequence the energy-separated (fiux unfolded) cross section obtained for the Mini- 
BooNE experiment while being higher than that obtained from SuSA alone still underestimates the 
experimental result even when 2p2h contributions are added. Recently, a strong difference between 
neutrino and antineutrino cross sections has been obtained within this model, with the 2p2h effects 
being significantly larger for antineutrinos than for neutrinos |52j . 

Two other effective models to account for MEC/2p2h effects have been proposed by Bodek et 
al. [S3] [transverse enhancement model (TEM)] and Lalakulich et al. [S3]. The TEM can easily be 
implemented in MC event generators [55]. It assumes that it is sufficient to describe properly an 
enhancement of the transverse electron QE response function keeping all other ingredients as in the 
free nucleon target case. Thus, some effective proton and neutron magnetic form factors are fitted to 
electron- nucleus data and later they are used, together with the free nucleon axial current, to study 
CCQE processes. It is to say, the TEM assumes that there are no nuclear medium effects (RPA, 2p2h 
mechanisms, etc..) affecting those nuclear response functions induced by the nucleon axial- vector 
current. Despite of a certain phenomenological success to describe the MiniBooNE data [531 155], such 
assumption seems quite unjustified. 

In the model of Ref. [54) . the multinucleon mechanism contributions are parametrized as phase 
space multiplied by a constant which is fitted to the difference of the energy-separated MiniBooNE 
data and the calculated QE cross section. RPA effects are not taken into account in [S3]. Since these 
tend to lower the cross section in particular at forward muon angles, the model of |^ underestimates 
the contributions of 2p2h effects there. Indeed, the authors of this reference find that the shape and 
over-all size of the 2p2h contribution turns out to be rather independent of the muon angle. This is in 
sharp contrast with the microscopical results obtained within the IFIC [iQ] [31] and SuSa models [52] , 
that find the 2p2h contribution becomes significantly less important as the muon scattering angle 
increases. 

Perspectives to measure the MEC/2p2h contribution 

The unambiguous experimental measurement of the MEC contribution to the CC inclusive cross sec- 
tion can be made by detecting hadrons in the final state. All the microscopic models provide up to now 
only the MEC/2p2h contribution to the muon inclusive 2D differential cross section: cPauie/dQ,{k')dE'^. 
Such models cannot describe detailed exclusive cross sections (looking into the nucleon side), as ex- 
plicit FSI effects, that modify the outgoing nucleon spectra, have not been addressed yet in these 
microscopical models. It is reasonable to assume that at the level of the primary reaction mechanism, 
they produce only slightly changes in cPai^^i/dQ{k')dE'^, leaving almost unchanged the integrated cross 
sections [22] [23]. 

A model to describe hadrons in the final state was proposed in |55J. It was implemented in 
the NuWro MC event generators and its predictions were used in the analysis of recent MlNERvA 
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antineutrino CCQE data. 

In the papers [Ml ES] various observables are discussed which can be used to detect MEC contribu- 
tion. One option is to look at proton pairs in the final state. Another possibility is to investigate the 
distribution of visible energy which allows to include contributions from protons below reconstruction 
threshold. The basic intuition from the electron scattering is that MEC events populate the region 
between QE and A peaks. Typically, to have a MEC event more energy must be transfered to the 
hadronic system than for a CCQE one. However, it should be stressed that the precision with which 
FSI effects are currently handled in MC codes can make such a measurement difficult. During last 
few years FSI studies were focused on pions only [56] aiming at understanding recent pion production 
data on nuclear targets |57j . Nucleons in the final state were never studied with a similar precision so 
there is less data to benchmark nucleon FSI effects. 



0.3.4 Monte Carlo event generators 

Monte Carlo codes (GENIE, NuWro, Neut, Nuance, etc) describe CCQE events using a simple RFC 
model, with FSI effects implemented by means of a semi- classical intranuclear cascade. NuWro offers 
also a possibility to run simulations with spectral function and an effective momentum dependent 
nuclear potential. It is also by now the only MC generator with implementation of MEC dynamics. 
Since the primary interaction and the final state effects are effectively decoupled, FSI do not change 
the total and outgoing lepton differential cross sections. 



0.3.5 Neutrino energy reconstruction 

Neutrino oscillation probabilities depend on the neutrino energy, unknown for broad fluxes and often 
estimated from the measured angle and energy of the outgoing charged lepton I only. This is the 
situation of the experiments with Cherenkov detectors where protons in the final state are usually 
below the Cherenkov threshold. Then, it is common to define a reconstructed neutrino energy Sj-ec 
(neglecting binding energy and the difference of proton and neutron masses) as: 

ME,-m\l2 

M - Ei + \pt\cosei ^ ' 

which would correspond to the energy of a neutrino that emits a lepton, of energy E(^ and three- 
momentum pi, with a gauge boson W being absorbed by a free nucleon of mass M at rest in a CCQE 
event. Each event contributing to the flux averaged double differential cross section da /dE^d cos Oi 
deflnes unambiguously a value of -Erec- The actual ("true") energy, E, of the neutrino that has 
produced the event will not be exactly E'rec- Actually, for each E^^c, there exists a distribution of true 
neutrino energies that give rise to events whose muon kinematics would lead to the given value of i^rec- 
In the case of genuine QE events, this distribution is sufficiently peaked (the Fermi motion broadens 
the peak and binding energy shifts it a little) around the true neutrino energy to make the algorithm in 
Eq. ([T]) accurate enough to study the neutrino oscillation phenomenon [58j or to extract neutrino flux 
unfolded CCQE cross sections from data (assuming that the neutrino flux spectrum is known) [591160]. 
The effect of this assumption on the much more demanding measurement of CP-violation effects is 
currently being evaluated. 

However, and due to presence of multinucleon events, there is a long tail in the distribution of true 
energies associated to each Ey-ec that makes the use of Eq. ([T]) unreliable. The effects of the inclusion 
of multinucleon processes on the energy reconstruction have been noticed in |55| and investigated in 
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Ref. |59j . within the Lyon 2p2h model and also estimated in Ref. [H], using the simplified model of 
Ref. [5l| for the multinucleon mechanisms. This issue has been more recently also addressed in the 
context of the IFIC 2p2h model in Ref. [60J , finding results in a qualitative agreement with those of 
Refs. p] and f6l]. 

In Ref. [60j it is also studied in detail the ^^C unfolded cross section published in [S]. It is 
shown there that the unfolding procedure is model dependent. Moreover, it is also shown that the 
MiniBooNE published CCQE cross section as a function of neutrino energy differs from the real one. 
This is because the MiniBooNE analysis assumes that all the events are QE. The authors of [60] finally 
conclude that the MiniBooNE unfolded cross section exhibits an excess (deficit) of low (high) energy 
neutrinos, which is mostly an artifact of the unfolding process that ignores multinucleon mechanisms. 

0.3.6 NC elastic 

MiniBooNE has also measured flux integrated NC elastic reaction cross-section [62]. Using these 
data, the best fit value of the axial mass was found te be Ma = 1.39 it 0.11 GeV. The measurement 
was possible because the MiniBooNE Cherenkov detector can observe also scintillation light from low 
momentum nucleons. An attempt was done to measure the nucleon strange quark component using 
the proton enriched sample of events with a result consistent with zero: As = 0.08 it 0.26. 

Theoretical considerations 

The MiniBooNE NCEl data were analyzed in [63j . The fit was done to the distribution of events 
with the best fit value of Ma equal to 1.28 it 0.05 GeV. Moreover the authors of concluded that 
axial mass as large as 1.6 GeV is still too small to reproduce the MiniBooNE NCEl data. Critical 
discusson of this statement can be found in Ref. |65j . 

0.4 The Resonance Region 

In the RES region the degrees of freedom are hadronic resonances, with the most important being 
the A(1232). Typical final states are those with a single pion. During the last five years several new 
pion production measurements have been performed. In all of them the targets were nuclei (most 
often carbon) and interpretation of the data in terms of the neutrino-nucleon cross section needed 
to account for nuclear effects, impossible to do in a model independent manner. Because of that 
it has become a standard that the published data include nuclear effects with most uncertain FSI. 
Perhaps not surprisingly, in several papers old deuterium ANL and BNL pion production data were 
re-analyzed aiming to better understand the pion production reaction on free nucleons. Theoretical 
models became more sophisticated and the major improvement was a development of well justified 
mechanisms for the non-resonant contribution in the A region. Some papers addressed the problem 
of higher resonances, a topic which will be investigated experimentally with future MINERvA results. 
On the other hand, there has been a lot of activity in the area of the coherent pion production and 
this subject will be discussed separately. 
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0.4.1 Experimental Results 

NCvrO 

Neutral current vr*^ production (NCvr'') is a background to appearance oscillation signal. One is 
interested in a vr'^ leaving the nucleus and recent experimental data are given in this format with 
all the FSl effects included. Signal events originate mostly from: NClvr'^ primary interaction with 
a vr*^ not being affected by FSl and NClvr"'' primary interaction with the tt^ being transformed into 
vr*^ in a charge exchange FSl reaction. An additional difficulty in interpreting the NCvr" production 
comes from a coherent (COH) contribution. In the case of MiniBooNE flux neutrino-carbon reactions 
(< Ejj >^ 1 GeV) it is estimated to account for ~ 20% of signal events j66j . 

Four recent measurements of NCvr^ production (K2K [HT], MiniBooNE neutrinos, MiniBooNE 
antineutrinos [68], SciBooNE [69j) are complementary. They use three different fluxes: (K2K, Fermilab 
Booster neutrinos and anti-neutrinos) and three targets: H2O (K2K), CH2 (MiniBooNE) and C^H^ 
(SciBooNE) . MiniBooNE presented the results in the form of absolutely normalized cross-section while 
K2K and SciBooNE reported only the ratios a{NClTT^) / a{CC). There is an important difference in 
what was actually measured: K2K and MiniBooNE present their results as measurements of final 
states with only one vr'^ and no other mesons. SciBooNE defines the signal as states with at least 
one vr^ in the final state so that a contamination from Itt'^Itt^, 2t:^ and > 27r (with > Ivr") final 
states is included and its fraction can be estimated to be 17% [5?]. Final results are presented as flux 
averaged distributions of events as a function of the vr*^ momentum, and in the case of MiniBooNE 
and SciBooNE also as a function of the tt*^ production angle. 

CC7r+ 

MiniBooNE measured CC Ivr'*' production cross sections, where the signal is deflned as exactly one 
Tr"*" in the final state with no other mesons |70j. A variety of flux integrated differential cross sections, 
often double differential were reported in and flnal state particles momenta. Also absolute vr"*" 
production cross sections as a function of neutrino energy are provided in Ref. [70]. The cross section 
results are much larger than NUANCE MC predictions and the difference is on average 23%. In Fig. 
[5] on the left GiBUU and NuWro predictions for CCvr"^ are compared to the MiniBooNE data. 

CCvrO 

MiniBooNE measured also CC Itt*^ production cross sections. As before, the signal is deflned as exactly 
one tt" in the flnal state [TTj. Various differential distributions are available. There is a dramatic 
discrepancy between the measured CC Ivr^ production cross section as a function of neutrino energy 
and NUANCE MC predictions in the region of lower energies. On average the data is larger by 
56 lb 20%, but for < 1 GeV the disagreement is as large as a factor of 2. In Fig. [s] on the right 
GiBUU predictions for CCtt^ are compared to the MiniBooNE data. 

Ratio a{CClTT+)/(j{CCQE) 

Another useful MiniBooNE measurement was the ratio a{CCl'K^) / a[CCQE) [72]. The ratio of 
CCl7r"'"-like (one pion in the flnal state) to CCQE-like cross-sections on CH2 as a function of neutrino 
energy was measured with an accuracy of ~ 10% in bins with highest statistics. This measurement 
puts constraints on the theoretical models which include QE, A excitation and MEC/2p2h dynamics. 
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Figure 5: Left: Differential cross section for CClTr"*" production in the final state (all the FSI effects are included). Mini- 
BooNE measurement [70] is compared to GiBUU [73] and NuWro computations. Right: the same for CCtt" production, 
but only GIBUU results are shown. 



But still, in order to compare with theoretical model predictions to these data, FSI effects must be 
included. In order to make such a comparison easier, MiniBooNE provided also FSI corrected data 
representing the ratio of CCl7r"'"/CCQE cross sections at the primary interaction. The corrected results 
are biased by MC assumptions and in particular they neglect most of the MEC/2p2h contributions 
which is contained in the QE-like sample of events. Finally, MiniBooNE re-scaled their results in order 
to get data points for an isoscalar target and enable comparison to old ANL and also more recent 
K2K data [73j . K2K measured ratio of cross sections on bound nucleons inside the nucleus corrected 
for FSI effects. CClvr"'" events were not identified on an event-by-event basis. 



Theoretical Considerations 

Due to nuclear effects a comparison to the new data is possible only for MC event generators, sophis- 
ticated computation tools like GiBUU and also a few theoretical groups which are able to evaluate 
FSI effects. 

Most of the interesting work was done within GiBUU. It turned out to be very difficult reproduce 
the MiniBooNE CClvr^ and CCItt*^ results: the measured cross section is much larger than theoretical 
computations. In the case of CC 1tt~^ production the discrepancy is as large as 100%. It was also 
noted that the reported shape of the distribution of vr"*" kinetic energies is different from theoretical 
calculations and does not show a strong decrease at > 120 MeV located in the region of maximal 
probability for pion absorption. 

The authors of [73| mention three possible reasons for the data/GiBUU predictions discrepancy: 
(i) the fact that A excitation axial form factor was chosen to agree with the ANL data only, neglecting 
the larger cross section BNL measurements; (ii) hypothetical 2p-2h-l7r pion production contribution 



analogous to 2p-2h discussed in the Sect. 0.3.3 (iii) flux underestimation in the MiniBooNE experi- 
ment. For the last point, the argument gets support from the better data/theory agreement found for 
the ratio, as discussed below. 

In the case of NCtt" production, a systematical comparison was done with NuWro MC predictions 
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with an updated FSl model for pions [57]. The overall agreement is satisfactory. Shapes of the 
distributions of final state vr'^'s are affected by an interplay between pion FSl such as absorption and 
formation time effects, understood here as an effect of a finite A life-time. It is argued that NCvr'^ 
production data can be very useful for benchmarking neutrino MC event generators. 

Because of the apparent data/MC normalization discrepancy for the CC vr"*" production, the inter- 
esting data is that for the ratio a{CClTr^ — like) / (t{CCQE — like). This observable is free from the 
overall fiux normalization uncertainty. However, it is not a direct observable quantity because in the 
experimental analysis it is necessary to reconstruct the neutrino energy and the procedures applied 
for the denominator and numerator are different. Three theoretical predictions for the ratio were 
published. The Giessen group compared to the MiniBooNE ratio data using the model described in 
[75] with the FSl effects modeled by the GiBUU code [76] . There is a significant discrepancy between 
the model and the data points: the calculated ratio is smaller. For the K2K data, the GiBUU model 
computations are consistent with the experimental results. 

The a{CCl'K'^) / (j{CCQE) ratio was also analyzed in Ref. [77]. In this analysis many nuclear 
effects were included: the in medium A self-energy (both real and imaginary parts), FSl effects within 
the cascade model of Ref. [78j, RPA corrections for the CCQE... Computations did not include 
contributons from the non-resonant background and from higher resonances. The contribution from 
the coherent pion production evaluated with the model of Ref. |79| (about 5% of the vr^ signal, 
a surprisingly large fraction) was also included in computations. The model predictions agree with 
MiniBooNE measurement for Ey <\ GeV and are below MiniBooNE data for larger neutrino energies. 

Finally, NuWro MC results for the ratio given in Ref. ^80j are slightly below the data points for 
larger neutrino energies. 

0.4.2 Theoretical Analyses 

It has been known since ANL and BNL pion production measurements that although being a dominant 
mechanism, A excitation alone cannot reproduce the data and that nonresonant background terms 
must be included in the theoretical models. There were many attempts in the past to develop suitable 
models but usually they were not very well justified from the theoretical point of view. 

Nonresonant background 

A general scheme to analyze weak pion production in the A region based on the chiral symmetry 
was proposed a few years ago in [81]. The model is supposed to work well in the kinematical region 
W < 1.3 — 1.4 GeV i.e. in the A region. The background contribution is particulary important at 
the pion production threshold, for values of W near M + m^. Vector form factors are taken from 
the electroproduction data and fits to helicity amplitudes [82]. Although particulary important for 
the channels t'^n — )• i~pTT^ and t'^n — )• i~mr~^ the background terms contribute also to the channel 
v^p — >• i~pTT^ changing the fitted values of the nucleon-A transition matrix elements. A comparison 
to existing NC pion production data was done as well and a good agreement was also found. An 
interesting question raised by the authors of [M] is that of unitarity. Their approach does not satisfy 
requirements of the Watson theorem and this can have some consequences e.g. worse agreement with 
the antineutrino pion production data. 

The model of the nonresonant background was used by the Giessen group which made several 
qualitative comparisons to both the ANL and BNL pion production data in the region W < lA GeV 
neglecting deuterium effects [HB]- In the case of neutron channels the model predictions are much 
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below the BNL data points and this is because the axial form factor parameters were optimized to 
the ANL data only. This choice goes back to the paper [82j where the authors came to the conclusion 
that the ANL and BNL data for the A"*""*" excitation are not compatible. 

Reanalysis of old bubble chamber data 

The issue of nucleon-A transition matrix element was discussed also in other papers. The questions 
are: what is the value of the 6*5^ (0)? How relevant are deuterium nuclear effects in dealing with ANL 
and BNL data? How much tension is there between both data samples? 

In Ref. [81] a fit was done to the ANL data in the A^^ channel only with the results: C^{0) = 
0.867 ± 0.075 and Maa = 0.985 ± 0.082 GeV. The obtained value of C^(0) was very different from 
what follows from off-diagonal Goldberger-Treiman relations {C^{0) ~ 1.15). 

The authors of [8D| made a fit to both ANL and BNL data including in the x^j terms with the 
overall flux normalization uncertainties, separate for ANL and BNL. In the fit the deuterium nuclear 
effects were included as correction factors to the distributions of events, using the results of |8l] . 
The main conclusion was that ANL and BNL data are in fact consistent. This statement was verified 
in a rigorous way using parameter goodness of fit method [85j • In the dipole parameterization of the 
C^(Q2) form factor the best fit values were found to be C^(0) = 1.14 ± 0.08 and Ma = 0.95 ± 0.04. 
Only A"'"^ channel was analyzed and like in Ref. |75j non-resonant background contributions was not 
included. 

So far the most complete analysis of both ANL and BNL data was performed in |86j : a nonresonant 
background was included and also deuterium effects were taken into account in the systematic way. 
The authors made several fits with various assumptions (see Table I) and in the fit IV they obtained 
C^(0) = 1.00 ±0.11. 

Other theoretical approaches 

In Ref. [87] the dynamical pion cloud effects are imposed on bare quark — A transition matrix 
elements. The model is able to reproduce both ANL and BNL weak pion production data. 

The authors of |88j focus on the consistent use of the A propagator. They show that the compu- 
tations relying on the standard Rarita-Schwinger propagator could lead to an underestimation of the 
weak pion production cross section. 

0.4.3 Coherent pion production 

In coherent pion production (COH) the target nucleus remains in the ground state. There are four 
possible channels, for CC and NC reactions, neutrinos and anti-neutrinos. A clear experimental signal 
for the COH reaction for high energies was observed and the aim of recent measurements was to 
fill a gap in the knowledge of a region around ~ 1 GeV COH cross-sections. At larger neutrino 
energies a recent measurement was done by MINOS which reported a NC reaction cross section at 
< Eu >= 4.9 GeV to be consistent with the predictions of the Berger-Sehgal model (see below). 

Experimental Results 

In the case of the NC reaction, MiniBooNE |66j and SciBooNE j89j searched for the COH component. 
SciBooNE [89] evaluated the ratio of the COH NCvr" production to the total CC cross-section as 
(1.16 ±0.24)%. 
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For the NC reaction MiniBooNE evaluated the COH component (plus possible hydrogen difFractive 
contribution about which little is known) in the NCvr'^ production as 19.5% (at < Ei, >~ 1 GeV) and 
then the overall flux averaged overall NClvr^ cross- section as (4.76 ± 0.05 ± 0.76) • 10~^°cmVnucleon. 
Unfortunately, it is difficult to translate both measurements into the absolutely normalized value of 
the NC COH cross-section because of strong dependence on the NUANCE MC generator used in the 
data analysis. In NUANCE, RES, COH and BGR (nonresonance background) NCvr*' reactions are 
defined according to primary interaction and COH pions are also subject to FSI. In the MiniBooNE 
analysis the fit is done for the composition of the sample of NCvr'^ events in terms of three components, 
and the COH fraction is defined as xcoH / {xcoH + xres)- 

In the case of the CC reaction, K2K [90] and SciBooNE [91] reported no evidence for the COH 
component. For the K2K analysis, the 90% confidence limit upper bound for the COH cross-sections 
on carbon was estimated to be 0.6% of the inclusive CC cross-section. The SciBooNE upper limits 
(also for the carbon target) are: 0.67% at < >~ 1.1 GeV, and 1.36% at < E^ >~ 2.2 GeV. 
SciBooNE reported also the measurement of the ratio of CC COH 7r+ to NC COH vr" production 
and estimated it as 0.14il^g |g. This is a surprisingly low value, which disagrees with results from 
the theoretical models which at SciBooNE energies typically predict values somehow smaller 2. For 
massless charged leptons isospin symmetry implies the value of 2 for this ratio and the finite mass 
corrections make the predicted ratio smaller. 

Theoretical developments 

Higher neutrino energy {E^ >~ 2 GeV) COH production data (including recent NOMAD mea- 
surement) were successfully explained with a PCAC based model [92]. Adler's theorem relates 
<ycoH{v + X 1/ + X + 7r°) at to (T(7r° + X 7r° + X). Subsequently, the model for 

the CC reaction, has been upgraded [93] to include lepton mass effects important for low Ey studies. 
The new model predicts the (Tcoh{'^^) / (^coHiT^^) ratio at E^ = 1 GeV to be 1.45 rather than 2. 
Another important improvement was to use a better model for dai^n + ^^C — )■ tt + ^'^C)/dt in the 
region of pion kinematical energy 100 MeV< T^r < 900 MeV. As a result, the predicted COH cross 
section from the model became reduced by a factor of 2-3 [94J. The PCAC based approach was also 
discussed in [95] and critically re-derived in Ref. [96] . 

At lower energies the microscopic A dominance models for the COH reaction [97J are believed to be 
more reliable. Within microscopic models there are still various approaches e.g due to differences in the 
treatment of the nonresonant background. The absolute normalization of the predicted cross-section 
depends on the adopted value of the N ^ A form factor C^{0) because acoH ~ (^5^(0))^- 

0.4.4 MC generators 

Almost all MC events generators rely on the old Rein-Sehgal resonance model for pion resonance 
production [98]. The model is based on the quark resonance model and includes contributions from 
18 resonances covering the region W < 2 GeV. The model is easily implementable in MC generators 
and it has only one set of vector and axial form factors. In the original model, the charged lepton is 
assumed to be massless and prescriptions to cope with this problem were proposed in Refs. [931 199]- 
It was also realized that the RS model can be improved in the A region by modyfying both vector 
and axial form factors using either old deuterium or new MiniBooNE pion production data |18| llOOj . 

As for coherent pion production, all the MCs use the Rein-Sehgal COH model [92]. The analysis 
of of MC event generators and theoretical models done in show that in the 1 — 2 GeV energy 
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region, the Rein Sehgal COH model predictions disagree significantly with all the recent theoretical 
computations and experimental results. 

A crucial element of MC is the FSI model. These are typically semiclassical intra-nuclear cascade 
models. The topic of FSI goes far beyond the scope of this review and we only note that the progress 
in understanding the experimental data requires more reliable FSI models. The existing models should 
be systematically benchmarked with electro and photoproduction data as it was done in the case of 
GIBUU. 

0.4.5 Duality 

Bridging the region between RES and DIS (where with a good approximation interactions occur on 
quarks) dynamics is a practical problem which must be resolved in all MC event generators. In 
MC event generators "DIS" is defined as "anything but QE and RES", what is usually expressed 
as a condition on the invariant hadronic mass of the type W > 1.6 GeV or 1^ > 2 GeV or so. 
Notice however that such a definition of "DIS" contains a contribution from the kinematical region 
< 1 GeV^ which is beyond the applicability of the genuine DIS formalism. RES/DIS transition 
region is not only a a matter of an arbitrary choice but is closely connected with the hypothesis of 
quark-hadron duality. 

Investigation of structure functions introduced in the formalism of the inclusive electron-nucleon 
scattering led Bloom and Gilman to the observation that the average over resonances is approximately 
equal to the leading twist contribution measured in the completely different DIS region. One can 
distinguish two aspects of duality: (i) resonant structure functions oscillate around a DIS scaling curve; 
(ii) the resonant structure functions for varying values of slide along the DIS curve evaluated at 
fixed Ql)is- 

In order to quantify the degree in which the duality is satisfied one defines the ratio of integrals 
over structure functions from RES and DIS: 

n{Q^, Qlro) = i^-" ' (8) 

it:d^m,Qiis) 

The integrals are in the Nachtmann variable f (x, Q^) = == and the integration region is 

; 1+^1+4x2 M2/Q2 ^ 

defined as Wmin < W < Wmax- Typically Wmin = M + m,r and Wmax = 1-6, ...,2.0 GeV. In the case 
of DIS, the value of Q^js is much larger and as a consequence the integral over ^ runs over a quite 
different region in W . 

Neutrino-nucleon scattering duality studies are theoretical in their nature because the precise data 
in the resonance region are still missing. The duality was studied in three papers: [57 } IIOH I102j . For 
neutrino interactions the duality can be satisfied only for the isospin average target. This is because 
the RES structure functions for proton are much larger than for neutron and in the case of DIS 
structure functions the situation is opposite. 

Theoretical studies were done with a model which contains resonances from the first and second 
resonance region but not the background contribution and with the Rein-Sehgal model which is com- 
monly used in MC event generators. If the resonance region is confined to W < 1.6 GeV the duality 
as defined in Eq. ([s]) is satisfied at the 75-80% level. If the resonance region is extended toW<2 GeV 
the value of the integral in Eq. ([s]) is only about 50%. These results are to some extent model depen- 
dent but a general tendency is that for larger W, DIS structure functions are much larger than the 
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resonance contribution, as clearly seen from Fig 3 in |101j and Fig. 7 in |102j . As shown in [1021 there 
is also a 5% uncertainty coming from an arbitrary choice of Q^jg- 

Two component duality hypothesis states that resonance contribution is dual to the valence quarks 
and the nonresonant background to the sea. Investigation done within the Rein Sehgal model with W < 
2 GeV revealed no signature of two component duality. Quark-hadron duality was also investigated 
in the case of neutrino nucleus interactions |103j . 

As a practical procedure for addressing this region, Bodek and Yang |104j have introduced and 
refined a model that is used by many contemporary neutrino event generators such as NEUGEN and 
its successor GENIE to bridge the kinematic region between the Delta and full DIS. The model has 
been developed for both neutrino- and electron-nucleon inelastic scattering cross sections using leading 
order parton distribution functions and introducing a new scaling variable they call 

Non-perturbative eff'ects that are prevalent in the kinematic region bridging the resonance and 
DIS regimes are described using the scaling variable, in combination with multiplicative K factors 
at low Q^. The model is successful in describing inelastic charged lepton-nucleon scattering, includ- 
ing resonance production, from high-to-low Q^. In particular, the model describes existing inelastic 
neutrino-nucleon scattering measurements. 

Their proposed scaling variable, is derived using energy momentum conservation and assump- 
tions about the initial/final quark mass and Pt- Parameters are built into the derivation of to 
account (on average) for the higher order QCD terms and dynamic higher twist that is covered by an 
enhanced target mass term. 

At the juncture with the DIS region, the Bodek- Yang model incorporates the GRV98 |105) LO 
parton distribution functions replacing the variable x with They introduce "K- factors" , different for 
sea and valence quarks, to multiply the PDFs so that they are correct at the low photo-production 
limit. A possible criticism of the model is the requirement of using the rather dated GRV98 parton 
distribution functions in the DIS region so the bridge to the lower W kinematic region is seamless. 

0.5 u-A Deep-inelastic Scattering: Introduction 

Although deep-inelastic scattering (DIS) is normally considered to be a topic for much higher energy 
neutrinos, wide-band beams such as the Fermilab NuMI and the planned LBNE beams do have real 
contributions from DIS that are particularly important in feed-down to the background that must 
be carefully considered. In addition, there are x-dependent nuclear effects that should be taken into 
account when comparing results from detectors with different nuclei and even when comparing results 
from "identical" near and far detectors when the neutrino spectra entering the near and far detectors 
are different. 

For this review, the definition of deep-inelastic scattering (DIS) is the kinematic based definition 
with W > 2.0 GeV and > 1.0 GeV. This is mostly out of the resonance production region and 
allows a fit to parton distribution functions. As said in Introduction, this is unfortunately not the 
definition used by several modern Monte Carlo generators that do not differentiate between simply 
"inelastic" interactions and deep- inelastic interactions calling everything beyond the Delta simply DIS. 
This is an unfortunate confusing use of nomenclature by the generators. 

In general, deep-inelastic scattering offers an opportunity to probe the partonic structure of the 
nucleon both in its free state and when the nucleon is bound in a nucleus. Description of the partonic 
structure can include parton distribution functions (PDFs) giving the longitudinal, transverse and spin 
distributions of quarks within the nucleon as well as, for example, the hadron formation zone giving 
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the time/length it takes for a struck quark to fully hadronize into a strong- interacting hadron. 

Neutrino scattering can play an important role in extraction of these fundamental parton distri- 
bution functions (PDFs) since only neutrinos via the weak-interaction can resolve the flavor of the 
nucleon's constituents: v interacts with d, s, u and c while the V interacts with u, c, d and s. The 
weak current's unique ability to "taste" only particular quark flavors significantly enhances the study 
of parton distribution functions. High-statistics measurement of the nucleon's partonic structure, 
using neutrinos, could complement studies with electromagnetic probes. 

In the pursuit of precision measurements of neutrino oscillation parameters, large data samples 
and dedicated effort to minimize systematic errors could allow neutrino experiments to independently 
isolate aU SIX the weak structure functions F'(^{x,Q'^), Ff^{x,Q^), F^^{x,Q^), F^^{x,Q^), 
xF^^ {xjQ"^) and xF^^ {x^Q"^) for the first time. Then, by taking differences and sums of these 
structure functions, specific parton distribution functions in a given (x, Q^) bin can in turn be better 
isolated. Extracting this full set of structure functions will rely on the y-variation of the structure 
function coefficients in the expression for the cross-section. In the helicity representation, for example: 



G' 



dxdQ^ 2ttx 



]^{F!^{x,Q^) + xFi^{x,Q^)) + 



{F!^{x, Q2) _ xFi^^x, Q2)) _ 2y^Fl(x, Q^)] . (9) 

where Fl is the longitudinal structure function representing the absorption of longitudinally polarized 
Intermediate Vector Boson. 

By analyzing the data as a function of (1 — y)^ in a given (x, Q^) bin for both u and v, all six structure 
functions could be extracted. 

Somewhat less demanding in statistics and control of systematics, the "average" structure functions 
F2{x, Q^) and xF^i^x, Q^) can be determined from fits to combinations of the neutrino and antineutrino 
differential cross sections and several assumptions. The sum of the u and V differential cross sections, 
yielding F2 then can be expressed as: 



dxdy dxdy vr L V ^ 2E 2 1 + Rt, ) 2/ 



1 + Rl 



(10) 



where Rl is equal to aL / ctt and now F2 is the average of F2 and F2 and the last term is proportional 
to the difference in xi*3 for neutrino and antineutrino probes, AxFa = xF^ — xFJ^ . In terms of the 
strange and charm parton distribution function s and c, at leading order, assuming symmetric s and 
c seas, this is Ax (s — c). 

The cross sections are also corrected for the excess of neutrons over protons in the target (for 
example the Fe correction is 5.67%) so that the presented structure functions are for an isoscalar 
target. A significant step in the determination of F2{x,Q'^) in this manner that affects the low-x 
values is the assumed AxF^ and Rl{x,Q'^). Recent analyses use, for example, a NLO QCD model as 
input (TRVFS |106) ) and assumes an input value of Rl{x,Q'^) that comes from a fit to the world's 
charged-lepton measurements jl07j . This could be an additional problem since, as will be suggested, 
Rl{x,Q'^) can be different for neutrino as opposed to charged-lepton scattering. 

The structure function xF^ can be determined in a similar manner by taking the difference in v 
and T' differential cross sections. 
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0.5.1 The Physics of Deep-inlastic Scattering 

There have been very few recent developments in the theory of deep-inelastic scattering. The theory 
has been well-established for years. The most recent developments in neutrino DIS scattering involve 
the experimental determination of parton distribution functions of nucleons within a nucleus, so- 
called nuclear parton distribution functions (nPDF). The more contemporary study of v nucleus 
deep-inelastic scattering using high-statistics experimental results with careful attention to multiple 
systematic errors began with the CDHSW, CCFR/NuTeV u Fe, the NOMAD u C and the CHORUS 
v Pb experiments. Whereas NuTeV |108j and CHORUS |109j Collaborations have published their full 
data sets, NOMAD [llOj has not yet done so. This short summary of DIS physics will concentrate on 
nuclear /nucleon parton distribution functions. 

Low-and-High Structure Functions: Longitudinal and Transverse 

Since the current and future neutrino beams designed for neutrino oscillation experiments will be 
concentrating on lower energy neutrinos (1-5 GeV), many of the interactions will be at the lower-Q 
edge of DIS or even in the "soft" DIS region - namely, W > 2.0 GeV however, with < 1.0 GeV^. 
Understanding the physics of this kinematic region is therefore important. 

Since both the vector and axial-vector part of the transverse structure function Ft go to at 
= (similar to £^ charged- lepton vector current scattering), the low-Q^ region v and 17 cross 
sections are dominated by the longitudinal structure function F^. The longitudinal structure function 
is composed of a vector and axial-vector component F^'-" and F^'-" and the low-Q^ behavior of these 
components is not the same as in the transverse case. The conservation of the vector current (CVC) 
suggests that Fj^'-' behaves as the vector current in charged-lepton scattering and vanishes at low 
Q^. However, the axial- vector current is not conserved and is related to the pion field via PCAC, 
so there is a surviving low contribution from this component [lllj and F^'-' dominates the low 
behavior. Consequently, the ratio R = Fl/Ft is divergent for neutrino interactions. This is 
substantially different from the scattering of charged leptons for which R is vanishing as and using 
measurement of R from charged lepton scattering to determine F2 for neutrino scattering is obviously 
wrong for lower Q. In addition, this non-vanishing and dominant longitudinal structure function could 
be important for the interpretation of low-Q^ nuclear effects with neutrinos to be described shortly. 

Low-and-High Structure Functions: l/Q^ Corrections 

Using a notation similar to that of reference |112j . the total structure function can be expressed in a 
phenomenological form: 



where i = 1, 2, 3 refers to the type of the structure function. Using i = 2 as an example, then F^^ is 
the leading twist component that has already included target mass corrections (TMC) and C4 is the 
coefficient of the twist-4 term, the first higher- twist term proportional to l/Q"^). There are, of course, 

further higher-twist terms ^ Qi + ■ ■ ■ proportional to ever increasing powers of l/Q however, for 
most phenomenological fits, the dominant leading twist plus twist-4 term are sufficient to describe the 
data. The target mass corrections are kinematic in origin and involve terms suppressed by powers 
of M^/Q^ while the higher twist terms are dynamical in origin and are suppressed as mentioned by 




(11) 
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powers of \/Q^. These higher- twist terms are associated with multi-quark or quark and gluon fields 
and it is difficult to evaluate their magnitude and shape from first principles. As with the kinematic 
target mass corrections, these must be taken into account in analyses of data at low and especially 
at large x. At higher the contribution of the HT terms is negligible and there are various global 
fits [1131 1114j to the structure functions (among various scattering input) to determine the parton 
distribution functions (PDFs) that do not include any HT terms. 

The analysis of nuclear PDFs to be described shortly uses data from a TeVatron neutrino experi- 
ment at very high neutrino energies and thus is one of the analyses that does not need to be concerned 
with higher-twist corrections. However, the current neutrino-oscillation oriented beam-lines are not 
high-energy and the analyses of this data may indeed need to consider both target mass corrections 
and higher- twist. If indeed inclusion of higher- twist in these analyses becomes necessary, the authors 
of [112j stress the importance of explicitly including both the target mass corrections and the higher 
twist corrections, even though they have very different physical origin and can have very different x 
dependence. It is important to note, as mentioned, that there are both nucleon and nuclear PDFs 
depending on the target. The relation between them, called nuclear correction factors, are currently 
being studied for both u-A and £^A. There are early indications that the nuclear correction factors 
for these two processes may not be the same. 

0.6 Recent DIS measurements: Neutrino Iron Scattering Results 

The difficulty, of course, is that modern neutrino oscillation experiments demand high statistics which 
means that the neutrinos need massive nuclear targets to acquire these statistics. This, in turn, 
complicates the extraction of free nucleon PDFs and demands nuclear correction factors that scale 
the results on a massive target to the corresponding result on a nucleon target. The results of the 
latest study of QCD using neutrino scattering comes from the NuTeV experiment |1U8| . The NuTeV 
experiment was a direct follow-up of the CCFR experiment using nearly the same detector as CCFR 
but with a different neutrino beam. The NuTeV experiment accumulated over 3 million and 17 
events in the energy range of 20 to 400 GeV off a manly Fe target. A comparison of the NuTeV 
results with those of CCFR and the predictions of the major PDF-fitting collaborations (CTEQ and 
MRST dm El] ) is shown in Figure [6} 

The main points are that the NuTeV F2 agrees with CCFR for values of xbj < 0.4 but is systemat- 
ically higher for larger values of xbj culminating at xbj = 0.65 where the NuTeV result is 20% higher 
than the CCFR result. NuTeV agrees with charged lepton data for xbj < 0.5 but there is increasing 
disagreement for higher values. Although NuTeV F2 and xF^ agree with theory for medium x, they 
find a different behavior at small x and are systematically higher than theory at high x. These 
results can be summarized in four main questions to ask subsequent neutrino experiments: 

• At high X, what is the behavior of the valence quarks as x — )• 1.0? 

• At all X and Q^, what is yet to be learned if we can measure all six u and u structure functions 
to yield maximal information on the parton distribution functions? 

• At all X, how do nuclear effects with incoming neutrinos differ from nuclear effects with incoming 
charged leptons? 

This last item highlights an overriding question when trying to get a global view of structure 
functions from both neutrino and charged-lepton scattering data. How do we compare data off nuclear 
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targets with data off nucleons and, the associated question, how do we scale nuclear target data to 
the comparable nucleon data. In most PDF analyses, the nuclear correction factors were taken from 
£='=-nucleus scattering and used for both charged-lepton and neutrino scattering. Recent studies by a 
CTEQ-Grenoble-Karlsruhe collaboration (called nCTEQ) jll6] have shown that there may indeed be 
a difference between the charged-lepton and neutrino correction factors. 

The data from the high-statistics Z/-DIS experiment, NuTeV summarized above, was used to per- 
form a dedicated PDF fit to neutrino-iron data |117j . The methodology for this fit is parallel to that 
of the previous global analysis [118| but with the difference that only Fe data has been used and no 
nuclear corrections have been applied to the analyzed data; hence, the resulting PDFs are for a proton 
in an iron nucleus - nuclear parton distribution function^ 

By comparing these iron PDFs with the free-proton PDFs (appropriately scaled) a neutrino-specific 
heavy target nuclear correction factor R can be obtained which should be applied to relate these two 
quantities. It is also of course possible to combine these fitted nPDFs to form the individual values of 
the average of F2{i'A) and F2{i'A) for a given x, to compare directly with the NuTeV published 
values of this quantity. This was recently done and the nCTEQ preliminary results |120j for low-Q^ 
are shown in Figure [7| Although the neutrino fit has general features in common with the charged- 
lepton parameterization, the magnitude of the effects and the x-region where they apply are quite 
different. The present results are noticeably flatter than the charged-lepton curves, especially at low- 
and moderate-x where the differences are significant. The comparison between the nCTEQ fit, that 
passes through the NuTeV measured points, and the charged-lepton fit is very different in the lowest-x, 
lowest-Q^ region and gradually approaches the charged-lepton fit with increasing Q^. However, the 
slope of the fit approaching the shadowing region from higher x where the NuTeV measured points and 
the nCTEQ fit are consistently below the charged-lepton A fit, make it difficult to reach the degree of 
shadowing evidenced in charged-lepton nucleus scattering at even higher Q^. 

The general trend is that the anti-shadowing region is shifted to smaller x values, and any turn- 
over at low X is minimal given the PDF uncertainties. More specifically, there is no indication of 
"shadowing" in the NuTeV neutrino results at low-Q^. In general, these plots suggest that the size of 
the nuclear corrections extracted from the NuTeV data are smaller than those obtained from charged 
lepton scattering. 

0.6.1 Comparison of the i'^A and uA Nuclear Correction Factors 

For the nCTEQ analysis, the contrast between the charged-lepton {(^A) case and the neutrino (I'A) 
case is striking. While the nCTEQ fit to charged-lepton and Drell-Yan data generally align with 
the other charged-lepton determinations, the neutrino results clearly yield different behavior as a 
function of x, particularly in the shadowing/anti-shadowing region. In the V case, these differences 
are smaller but persist in the low-x shadowing region. The nCTEQ collaboration emphasize that both 
the charged-lepton and neutrino results come directly from global fits to the data, there is no model 
involved. They further suggest that this difference between the results in charged-lepton and neutrino 
DIS is reflective of the long-standing "tension" between the light-target charged-lepton data and the 
heavy-target neutrino data in the historical global PDF fits |12H I122j . Their latest results suggest 
that the tension is not only between charged-lepton light-target data and neutrino heavy-target data, 
but also between neutrino and charged-lepton heavy-target data. In other words a difference between 

*For more details of the fitting tecliniques and resulting comparisons with charged-lepton scattering see Part II of 
reference [119| . 
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charged-lepton {i A) and the neutrino (I'A) when comparing the same A. 

Concentrating on this interesting difference found by the nCTEQ group, if the nuclear corrections 
for the £^A and uA processes are indeed different there are several far-reaching consequences. Con- 
sidering this, the nCTEQ group has performed a unified global analysis |116j of the £^A, DY, and 
vA data (accounting for appropriate systematic and statistical errors) to determine if it is possible 
to obtain a "compromise" solution including both £^A and uA data. Using a hypothesis-testing cri- 
terion based on the distribution that can be applied to both the total as well as to the of 
individual data sets, they found it was not possible to accommodate the data from uA and £^A DIS 
by an acceptable combined fit. 

That is, when investigating the results in detail, the tension between the £^Fe and uFe data sets 
permits no possible compromise fit which adequately describes the neutrino DIS data along with the 
charged-lepton data and, consequently, £^Fe and uFe based on the NuTeV results, have different 
nuclear correction factors. 

A compromise solution between vA and £^A data can be found only if the full correlated systematic 
errors of the vA data are not used and the statistical and all systematic errors are combined in 
quadrature thereby neglecting the information contained in the correlation matrix. In other words 
the larger errors resulting from combining statistical and all systematic errors in quadrature reduces 
the discriminatory power of the fit such that the difference between uA and £^A data are no longer 
evident. This conclusion underscores the fundamental difference [116j of the nCTEQ analysis with 
other contemporary analyses. 

On the other hand, a difference between uA and £^A is not completely unexpected, particularly 
in the shadowing region, and has previously been discussed in the literature \123\ I124j . The charged- 
lepton processes occur (dominantly) via 7-exchange, while the neutrino-nucleon processes occur via 
1^^-exchange. The different nuclear corrections could simply be a consequence of the differing prop- 
agation of the hadronic fluctuations of the intermediate bosons (photon, W) through dense nuclear 
matter. Furthermore, since the structure functions in neutrino DIS and charged lepton DIS are dis- 
tinct observables with different parton model expressions, it is clear that the nuclear correction factors 
will not be exactly the same. What is, however, unexpected is the degree to which the R factors differ 
between the structure functions F^'^^ and i^l^*^. In particular the lack of evidence for shadowing in 
neutrino scattering at low down to x ~ 0.02 is quite surprising. 

Should subsequent experimental results confirm the rather substantial difference between charged- 
lepton and neutrino scattering in the shadowing region at low-Q^ it is interesting to speculate on the 
possible cause of the difference. A recent study of EMC, BCDMS and NMC data by a Hampton 
University - Jefferson Laboratory collaboration |125j suggests that anti-shadowing in charged-lepton 
nucleus scattering may be dominated by the longitudinal structure function F^. As a by-product of 
this study, their figures hint that shadowing in the data of EMC, BCDMS and NMC fi A scattering 
was being led by the transverse cross section with the longitudinal component crossing over into the 
shadowing region at lower x compared to the transverse. 

As summarized earlier, in the low-Q^ region, the neutrino cross section is dominated by the longi- 
tudinal structure function F^ via axial-current interactions since Ft vanishes as as — t- similar 
to the behavior of charged lepton scattering. If the results of the NuTeV analysis are verified, one 
contribution to the different behavior of shadowing at low-Q^ demonstrated by A and £ A, in ad- 
dition to the different hadronic fiuctuations in the two interactions, could be due to the different mix 
of longitudinal and transverse contributions to the cross section of the two processes in this kinematic 
region. 



26 



Neutrino Physics 



Neutrino-Nucleus Interactions 



Another hypothesis of what is causing the difference between neutrino and charged-lepton shad- 
owing results comes from V. Guzey [126j who speculates that at low x, low-Q^ the value of y is 
close to unity and the neutrino interactions primarily probe the down and strange quarks. This is 
very different than the situation with charged-lepton scattering where the contribution from down and 
strange quarks are suppressed by a factor of 1/4 compared to the up and charm. Therefore, the dis- 
crepancy between the observed nuclear shadowing in lepton-nucleus total cross section at small x and 
shadowing in total neutrino-nucleus cross section could be caused by the absence of nuclear shadowing 
in the nuclear strange quark distributions as extracted from the neutrino-nucleus data or even the 
poor knowledge of the strange-quark distribution in the free-nucleon that affects the neutrino-nucleus 
ratio more than the charged-lepton ratio. 
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Figure 6: A comparison of the measurements of the F2 structure function by NuTeV and CCFR and 
the predictions from the global PDF fits of the MRST and CTEQ collaboration |115j that does not 
use the NuTeV data points as input to their fit. The model predictions have already been corrected 
for target mass and, most significantly, nuclear effects assuming these corrections are the same for 
charge-lepton and neutrino interactions 
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Figure 7: Nuclear correction factor R for the average F2 structure function in charged current I'Fe 
scattering at =1.2, 2.0, 3.2 and 5.0 GeV^ compared to the measured NuTeV points. The green 
dashed curve curve shows the result of the nCTEQ analysis of A (CHORUS, CCFR and NuTeV) 
differential cross sections plotted in terms of the average -P/''^ divided by the results obtained with the 
reference fit (free-proton) PDFs. For comparison, the nCTEQ fit to the charged-lepton data is shown 
by the solid blue curve. 
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